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1 Introduction 



Over the past few years there has been much work done^ 1 ^ studying the radiative energy loss of 
high energy partons passing through hot and cold matter. These studies are extensions to QCD 
of the analogous QED problem considered long ago by Landau, Pomeranchuk and Migdal B eI . In 
one version of this work initiated in Refs. Jl| and p] and continued in Refs. || and Q one follows 
the multiple scattering i'0 of the high energy parton in the QCD matter and the radiative 
gluon spectrum induced by the multiple scattering is evaluated. A number of interesting, and 
perhaps surprising results were found. The energy loss of a high energy jet in a hot QCD plasma 
appears to be much larger than that in cold nuclear matter even at moderate, say 200 MeV, 
temperatures of the plasma. When a very high energy parton passes through a length L of 
hot or cold matter the induced radiative energy loss is proportional to L 2 . A curious relation, 
—dE/dz = const • a s N c p 2 ±w , was found between the energy loss and the width of the Gaussian 
transverse momentum broadening of a high energy parton in QCD matter. 

A different and Yery elegant approach to the energy loss problem has been developed by 
Bronislav Zakharov oB. In his approach the gluon radiation probability is given by the difference 
of the probabilities of a virtual fluctuation of a high energy quark into a quark and a gluon 
occuring in the vacuum and in the medium. The formalism uses a clever device of treating 
the quark in the complex conjugate amplitude as an antiquark in the amplitude. What is 
finally calculated then is the amplitude for a quark-antiquark-gluon system to pass through a 
QCD medium without inducing inelastic reactions in the medium. For a medium having many 
scatterers such an amplitude is not small only if the quark-antiquark-gluon system is compact in 
transverse coordinate-space, and this compactness makes the process perturbatively calculable 
in QCD. 

One of the major purposes of the present work is to show that the BDMPS 'S § and Z i' 1 
formalisms are equivalent. However, before showing this equivalence it is necessary to extend the 
BDMPS formalism beyond the soft gluon approximation where it was originally formulated. It is 
also necessary to include some virtual graphs which were omitted in our original formulation and 
which led to some numerical discrepancies when our results were compared to those in Refs. |J 
and @. 

In Sees. 2 and 3 we outline a derivation of the BDMPS formulas for the radiative gluon 
spectrum for a high energy quark passing through either a hot or cold QCD medium, which 
formulas are valid even when the radiated gluon carries a finite fraction of the quark's momentum. 
Special care is taken to show which factors have changed from our previous results when all terms 
corresponding to virtual corrections (elastic scattering in the medium) are included. 

In Sec. 4 we evaluate our formulas for the radiative gluon energy spectrum. We consider both 
the case when the quark approaches the medium from outside and when the quark is produced 
by a hard scattering in a medium. 

In Sec. 5 we show the equivalence between our approach and that of B. Zakharov. We 
show that our formula (33) for the radiative spectrum leads to fl59|) which is equivalent to (4) of 
Ref. 0. We then explain more qualitatively why the two formalisms, seemingly very different, 
are in fact completely equivalent. 
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2 The Born term for the amplitude 



In this section we calculate the lowest order terms for the emission of a gluon from a quark which 
may enter the QCD matter from the vacuum or which may be produced in the matter through 
a hard interaction. We begin by giving the basic vertex for emission of a gluon of momentum k 
from a quark of momentum p. 

We shall continue to call the fermion from which the gluon is emitted a quark. However, in 
our final formulas results for spin 1/2 fermions of any colour representation, R, can be obtained 
by multiplying the gluon emission spectrum ( |42b|) by the ratio of the "colour charges", Cr/Cf- 
The generalisation to spin-1 (gluon) and spin-0 projectiles will also be given in Sec. |5\2 . 



2.1 The basic vertex 

The basic vertex for gluon emission is given by 



p-k, r 




u r i{p — k) 7 • e u r (p) 



x\/l — x 

We use a notation { } such that for any two-dimensional vector Vi, i = 1,2 

x\ x 
2 J 2 J 3 



e ■ { k — xp) . 



(2) 



where is the antisymmetric tensor; e\2 = —621 = 1, en = €22 = 0. The gluon polarisation 
vector is 

e • k e ■ k 



2k ' 2k 



In ([T|-0) we assume that 



(3) 

(4a) 
(4b) 

but we do not assume that x = k/p is small. The vertex is diagonal in the quark spin indices 
r and r' because of the vector nature of the coupling and because of our assumption that the 
fermion is massless. 



\k\ 


< 


k 


\p\ 


< 


p\ 



2.2 The Born terms for the amplitude; scattering in the medium 



In our procedure one explicitly integrates the emission time, t, of the gluon between the times 
the quark, or the quark-gluon system, scatters in the medium. Consider, for example, the graphs 
in Fig 2a where the quark scatters inelastically in the medium at time t\ and then later emits 
the gluon at time t. The t-dependent phase factor associated with this graph is 



exp 



it 



(p — k) 2 p 2 
2k ' 2(p - k) ~ 2p 



+ 



exp 



it 



(k — xp) 2 
2x(l — x)p 



(5) 
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The lower limit of the t-integration is t\ while the upper limit will be the time of the next inter- 
action or t = oo in case the quark-gluon system leaves the medium without further interactions. 
Keeping the lower limit of the integral of the expression (||) and multiplying by the gluon emission 
vertex gives the factor 



Aip\/l — x exp 



(k — xp) 2 
2x(l — x)p 



T a -e, 



(6) 



where 



{k — xp} 
(k — xp) 2 



(7a) 



Similarly one can evaluate the emission term and energy denominators in the remaining graphs 
b-g of Fig. 1 to obtain 



T 
—g 



{k — xp—(l— x)q] 
(k — xp — (l — x)q) 2 ' 
{ k — xp + xq} 



(k — xp + xq)' 



-T t 



{k — xp — q} 
(k — xp — q) 2 



(7b) 

(7c) 
(7d) 
(7e) 



The (—1) factors in some terms in fl7|) occur when t\ is the upper limit of the t integration, as 
happens for the graphs b, c, e, f and g of Fig. 1. 

Graphs a-c of Fig. 1 correspond to inelastic reactions with the medium while graphs d-g 
correspond to forward elastic scatterings in the medium. For terms a-c there are corresponding 
inelastic reactions in the complex conjugate gluon emission amplitude. In the approximation 
that the forward elastic amplitude for quark scattering off particles in the medium is purely 
imaginary, the elastic and inelastic terms are proportional to the same function V to be defined 
more precisely below. 

It is remarkable that the gluon emission vertices are proportional to exactly the same combi- 
nations of transverse momenta whose squares enter in the energy denominators for each of the 
terms in (|7]). This is a consequence of gauge invariance. The fact that the emission amplitude is 
inversely proportional to the first_power of the relevant transverse momentum follows from the 
Gribov bremsstrahlung theorem 0, the generalisation of the Low-Barnett-Kroll theorem @ to 
hard collinear radiation, x ~ 1, k±_ <C p. This will prove crucial in obtaining a simple form for 
the gluon emission spectrum when x is not small. 

It is our convention that the time of emission of the gluon in the complex conjugate amplitude, 
t f , is later than the emission time in the amplitude. The opposite sequence of times will be 
accounted for by a factor of 2 in the gluon spectrum (|31[). This convention means that for the 
terms in the complex conjugate amplitude corresponding to a-c of Fig. 1 it is only the quark 
and not the quark-gluon system which scatters inelastically off the medium at time t\. Thus the 
colour factors associated with the scattering in the medium for the graph in Fig. la is obtained 
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by considering the graph of Fig. 2a where the part of the graph to the right of the vertical line 
(cut) is the complex conjugate amplitude. The initial and final p-lines have the same colour. 
The graph in Fig. 2b corresponds to the graph in Fig. lb along with the corresponding complex 
conjugate amplitude. Using the formulas for colour factors given in Appendix B of Ref. || it 
is straightforward to evaluate the "colour factors" associated with the hooking of the g-lines in 
Fig. 1 (and Fig. 2) with the quark or quark-gluon passing through the medium, namely 



F a — — 2F d — — 2Ff — C F 

2 ' 



F b = -F e = F g = , (8b) 



F c = ——. (8c) 

In addition to pure colour factors the F's also include a factor of (—1) for virtual terms and a 
factor of 1/2 for virtual terms where the q and — q lines attach to the same quark or gluon line. 
These factors naturally appear in a Feynman diagram description of Glauber multiple scattering. 
For quark propagating through a QCD medium these factors guarantee probability conservation. 
Indeed, the virtual corrections to the amplitude and the complex conjugate amplitude, 2 * (— |), 
cancel the +1 coming from the scattering in the medium. 
It is convenient to introduce scaled momentum variables 

k/fj> = u, q/p = Q, p/p = Y_, (9) 

where /i is an appropriate scale for the problem. In case the high energy fermion moves through a 
hot QCD plasma \i may be taken to be the inverse Debye screening length, while for cold nuclear 
matter /i is naturally taken to be a typical transverse momentum exchanged in a quark-nucleon 
scattering. In addition to the factors and Fi the graphs in Fig. 1 naturally include a quark- 
"particle" cross section when corresponding complex conjugate amplitudes are included. Thus, 
for example, the graph of Fig. la when multiplied by a similar complex conjugate amplitude, 
leads to the graph illustrated in Fig. 2a, and this graph is clearly proportional to the quark- 
particle scattering cross section in a two-gluon-exchange approximation. We define a normalised 
quark-particle cross section o by 

V(Q 2 ) = — (10) 

7r a dQ z 

with 

° = j^Q- (ii) 

In case the medium is cold nuclear matter the "particle" referred to above is a nucleon while in 
case the medium is a hot QCD plasma the particle can be taken to be a quark or gluon. 

In this paper we do not consider collisional energy loss. To ensure the validity of the indepen- 
dent scattering picture and to guarantee that da/dQ 2 depend only on transverse momentum, it 
suffices to assume that the energy transfer from the quark to a particle in the medium be small 
compared to the incident energy II. 
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It is useful to combine factors to define a Born-term gluon emission amplitude fo(U_,V) as 



2C r 



F A 1 J 



(12) 



Here fo is an amplitude for gluon emission although it also embodies colour factors and quark- 
particle scattering factors from the complex conjugate amplitude. The factor N c /2Cf is included 
to agree with our earlier choice of normalisation B of fo(U_,V). 

We note that the Born amplitude fo and the full amplitude / to be introduced later in Sec. 0, 
as well as their impact-parameter images, fo{M) an d f(B_), are two-dimensional vectors. It is 
implied hereafter though we chose not to underscore /'s as transverse vectors. 

It is convenient to change from momentum space to impact parameter space. Since fo depends 
on the two momenta, U_ and V_, it might be expected that two impact parameters would be 
necessary. However, because k and p enter in (^|) only in the combination k — xp, the amplitude 
fo can depend on U_ and V_ only in the combination U_ — xV_. Thus it is possible to express fo in 
terms of a single "impact parameter" as 



fo(U,Y) 



d 2 B 
(2vr) 2 



JB-(U-xV) 



km 



(13) 



with 



km 



d 2 (U-xV) e 



-iB-(U_-xV_) 



fa(U,V). 



(14) 



Using 



IP 



2tt / B 2 



(15) 



it is straightforward to find 



km 



-2iri 



{§} 

B 2 



V(B(l-x)) 



V{B) 



1 



1 - V(-Bx) 



(16) 



with 



V(B) = j d 2 Q V(Q 2 ) ; V(0) = 1 . 



The { } symbol in ( p^6| ) is defined in (|2|). 

Now let us compare ( |16|) with our previous small-x results in (4.5) and (4.25) of Ref. 0. Of 
the three terms on the right hand side of (H^) the first and third terms come from the production 
terms in the medium, graphs a-c of Fig. 1, while the second term comes from the elastic scattering 
terms in the medium, graphs d-g of Fig. 1. In the small-x limit the first and second terms on 
the right hand side of (0) become equal while the third term vanishes. Comparing with (4.25) 
of Ref. [U we note that our present result is larger than what we previously found by a factor of 
2 in the small-x limit. This factor of 2 is due to an incomplete treatment of virtual contributions 
in our earlier work. 
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2.3 The Born amplitude for the hard scattering case 



In case the high energy quark is produced in the medium through a hard scattering one must 
also include the contribution coming when the endpoint of the integration of the gluon emission 
coincides with the time of the hard scattering. It is not important to know the details of the 
hard scattering since we are here only interested in the radiative gluon spectrum accompanying 
the hard scattering and subsequent rescatterings in the medium. Thus we may imagine that the 
gluon is produced in a collision of a highly virtual photon with a quark in the medium, which 
collision transfers a large energy to the struck quark which then rescatters in the medium and 
radiates a gluon. 




k 



The gluon emission amplitude is as given by T_ a in (ffa|). Since the gluon is emitted after the hard 
scattering there is no overall colour factor so that the basic Born term here is 



{U - xV} 

{u-xvy ' 



(17a) 



or 



- (17b) 



in impact parameter space. The expression given in ( |17b| ) must be added to t^/oGS), the 
Fourier transform of the left hand side of 



3 The time evolution of gluons in the medium 

After the gluon is emitted from the high energy quark, the quark-gluon system moves through 
the medium and carries out multiple scatterings i'0 with the particles of the medium. It is 
not certain that the gluon will be produced as a physical gluon until there is gluon emission in 
the complex conjugate amplitude. Thus, we must follow the time evolution of the quark-gluon 
system in the amplitude up to the time of emission in the complex conjugate amplitude. This is 
what will be done in this section. 

Let f(y_,Y_,t) be the quark-gluon amplitude at time t starting from fo(U_,V_) at t — 0. U_is 
the scaled gluon momentum while V_ — U is the scaled quark momentum as given in (||) for the 
graphs of Fig. 1. The amplitude f(U_,V_,t) is illustrated below. 
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We recall that t is less than the time t' at which the gluon is emitted in the complex conjugate 
amplitude. Then the time-dependence of / comes partly from the free evolution of the quark- 
gluon system and partly from interactions in the medium. Suppose t\ is the time of the last 
interaction with the medium before time t. Although / represents the gluon amplitude it also 
includes interactions of the quark with the medium in the complex conjugate amplitude as already 
noted for / . Then the possible interactions at t\ are shown in Fig. 3. For example the graph 
shown in Fig. 3a illustrates an inelastic interaction of the gluon with a particle in the medium in 
the amplitude and an interaction of the quark with the same particle in the complex conjugate 
amplitude. Perhaps the only unusual graph in Fig. 3 is graph e where the only interaction is a 
forward elastic scattering of the quark with a particle of the medium in the complex conjugate 
amplitude. 

The amplitude f(U_,Y_,t) obeys the integral evolution equation 



f(U,V,t) = e i4 ' it ~ tl) fo(U,V)+ [ dteW-«)£!L [ d 2 QV{Q 2 

Jh J 



N 

-jf{U-QjV-QjH) 



1 -mV-Q,t) - ^f{U,V,t') - f) - ^f{U,V,t') + ^f(U-Q,V,t') 



'11 



The first term in the right hand side of ( PD gives the free evolution of the quark-gluon system 
between t% and t in the amplitude and of the quark in the complex conjugate amplitude according 
to the phase factor 

e-*-> with , = £ + <^-g = (E-'K>V (19) 
Y 2k 2(p - k) 2p 2x(l - x)p 



Differentiating over t we obtain the equation 



N 

-£m-Q,v-Q,t) 



1 -f(U,V-Q,t) - ^f(U,V,t) - ^f(U,V,t) - ^f(U,V,t) + ^f(U-Q,V,t) 



2N C JK —'— — ' 2 Jy — '— ' ' 2 JX — '— ' ' 2 JV — '— ' ' 2 

(20) 

The first term in the right hand side of (^) comes from free propagation and the rest of the 
terms come from the upper limit of the integration over t\ at t\ = t. In (|T8D, (|20| ) p is the density 
of scatterers in the medium while a is the cross section of a quark with a particle in the medium. 
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The factor V(Q 2 )a corresponds to the differential cross section for scattering the high energy 
quark with momentum transfer Qfi while the 1 / Cp factor takes out the quark colour factor for 
quark scattering with a particle of the medium. The correct colour factors are then inserted in 
the various six terms in the integrand in fl2"0|). The six terms in the integrand in ( p0[ ) correspond, 
in order, to the six graphs of Fig. 3. The colour factors are 



F n = -F c 
F d 



p — _£f 



p . 1_ 

b ~ 2N C ' 



(21a) 
(21b) 
(21c) 



As in (H) we include (—1) factors in F for all virtual terms while factors 1/2 are included for the 
terms with both gluon lines, q and — q, hooking to the same high energy quark or gluon line. 

It is easy to see that f(U_,V_,t) = f(U_— x Y-it) since the explicit factors of U and V_ in (|20|) 
occur in the combination U_ — xV_ while f(U_,V_, 0) = fo(LL~ x Y.)- Introducing ®~ 



t N r 



A 2C f 



K 



A/x 2 



2C f 



A/x 2 



2C f 



2x(l - x)p N c 2u(l - x) N c 



(22) 



with pa — 1/A, where A is the mean free path of the high energy quark in the medium, one can 
write (|20f) as 

d 



dr 



f(U - xV, r) = ik(U - xV) 2 f(U - xV, r) 



- J d 2 Q V(Q 2 ) ( [f(U - xV, t) - f(U -xV-{\- x)Q, r)] (23) 
+ [f(U - xV, r) - f(U -xV-Q, t)] - [f(U - xV, r) - f(U -xV + xQ, r)] ^ . 
Now it is a simple matter to go to impact parameter space, 

f(B,r) = J <P(U-xV)e^ & '^®f(U-xV,T), 



(24) 



to find 



l-V(B(l-x)) 



V{B) 



1 

N 2 



l-V(-Bx) )f(B,r). 



(25) 



This is the basic Schrodinger-type evolution equation for the propagation of the quark-gluon 
system in a QCD medium. It should be solved with the initial condition 



imo) = urn 



(26) 



with fo(B_) given in (JT 



Comparing (|25|) to (4.23) of Ref. we again find a factor 2 discrepancy in the small-x limit 
of the second term on the right hand side of (^) as compared to (4.23). Our previous calculation 
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effectively amounted to keeping graphs a, b, d and e of Fig. 3. (Although graphs d and e were 
not explicitly considered, they were effectively included through the use of a mean free path 
term, the (—1) in the first term on the right hand side of (4.16) of Ref. ||.) That calculation 
was in error because the virtual terms, graphs c, d, e, f of Fig. 3 are not completely taken into 
account in the mean free path treatment of Ref. ||. We note that the coefficient of f{B_,r), in 
the second line of (^), the potential term, is of the same form as the three-body cross section 
used by Zakharov in eq. 23 of Ref. || . 



4 The Born term for the complex conjugate amplitude 

Now that we have calculated gluon emission in the amplitude and followed its evolution in time, 
it becomes necessary to calculate gluon emission in the complex conjugate amplitude. Suppose 
the gluon is emitted at time t in the complex conjugate amplitude. As usual we integrate this 
emission time between elastic or inelastic interactions with the medium. Except for colour factors 
the calculation proceeds exactly as in Sec. 2. Suppose t\, the time of interaction in the complex 
conjugate amplitude with the particle in the medium, serves as the endpoint of the t-integration. 
At t\ the amplitude consists of a quark-gluon system described by f{t\). The graphs describing 
the emission are shown in Fig. 4 where the vertical line indicates that the gluon is put on-shell. 
Terms to the left of the vertical line belong to the amplitude while to the right belong to the 
complex conjugate amplitude. We have rearranged the momenta so that the gluon emission 
amplitude f(U_ — xV_, t\) appears uniformly in all the graphs. 

The "colour factors" for gluon emission in the complex conjugate amplitude are 

\F a = -F b = F c = -F d = -F h = Fi = f , (27a) 
F e = -2F g = -2Fj = C F , (27b) 
Ff = -afc , (27c) 

where again a factor of (—1) for virtual terms and a factor of 1/2 for identical particles have 
been included in (0). The gluon emission terms, analogous to (0), for the complex conjugate 
amplitude are 



J ITT) r 

£ = Tl=r g = Tl = —Tj = - 1 -_ x - ) J 2 , (28a) 



( k — xp + xq) 
It is now straightforward to check that 



\k — xp — q}* 

s - = ( 28b ) 

~ c (k — xp—(l—x)q) 2 ' 

\k — xp + xq}* 

T* = 1= = (28d) 



i J r at 
— Y, J d 2 QV(Q 2 )T*(U,V,Q) Fi = -^-f*(U- x V), (29) 
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which, except for a minus sign, is similar to (0). Taking the Fourier transform gives 



J d\U-xV)e-W-*Z ! _Lj2 J d 2 QV(Q 2 )T*(U ,V,Q) F i = ft(B) ^j- , (30) 



with fo exactly the same function as given in (16) 



Comparing with what was found in Ref. Q, it is easy to verify that graphs a+b+c of Fig. 4 
generate what was called the U Y" term and that this contribution is equal to the first and second 
of the three terms on the right hand side of (jlB]). Graphs d+e+f generate what was called the 
U H" term which is equal to the second and third terms on the right hand side of flTBJ). The 
virtual graphs, graphs g+h+i+j of Fig. 4 were not included in Ref. || and it is easy to check 
using ( p7|) and (|28| ) that they give (—1) times the second term on the right hand side of (|i~6|) . 
Thus in the small- a; limit where the third term on the right hand side of (|16D is small we now 
find a factor 2, instead of the 3 coming from the sum of the Y + H graphs in Ref. ||, due to 
inclusion of the virtual terms. 



5 The spectrum of radiated gluons and energy loss 

In this section we remind the reader of the formula for the induced spectrum of radiated gluons fl . 
We then solve for the spectrum in the soft gluon limit for a volume of matter large enough 
that a high energy quark carries out many scatterings. Finally we integrate the spectrum of 
radiated gluons to find the radiative energy loss in circumstances where the energy loss problem 
is dominated by soft gluons, that is when a gluon having coherent length on the order of the 
dimensions of the medium has a longitudinal momentum much less that that of the high energy 
quark. 



5.1 The formula for the induced gluon radiation 

We give the formula for the induced gluon spectrum coming from a high energy quark produced 
in a hard collision in the medium. Then the Born term in the amplitude will be associated either 
with a scattering in the medium, as in (|1^), or with the hard vertex, as in ( |17a| ). In case the 
gluon is radiated from a high energy "on-shell" quark entering the medium one simply drops the 
term associated with the hard vertex. The induced gluon spectrum is 



du dz 7r 2 L 



c r ( r r iv jv 

^2Rej <i 2 {/|/ dh J dhpa^-fiU-xV^t-tJ-pa^-fZiU-xV) 



R-0 ( 31 ) 



The various terms in (|3~lD have simple interpretations. The a s Cp/7r 2 is the coupling of a gluon 
to a quark. The 1/L comes because we calculate the spectrum per unit length of the medium. 
The factor 

iV 

c ■ f{H- xV,t 2 -t 1 ) padt x 



2C f 



11 



gives the number of scatterers in the medium, pa dt\, times the gluon emission amplitude at ti, 
then evolved to t 2 . The factor 

-i^rfo(U-xV) padt 2 

gives the number of scatterers times gluon emission in the complex conjugate amplitude. The 
overall normalisation in the small-x limit is fixed by comparing with Ref. ||. We note that 
this normalisation is correct even when x is not small because the \/T — x present in the basic 
amplitude © but not included in our definition of fo cancels with a 1/1 — x quark phase space 
factor. 



The second term on the right hand side of (fHf ) gives the contribution of gluon emission due 



to the hard scattering in the medium. If there is no hard scattering in the medium this term 
should not be included. We assume the hard scattering happens at t — 0. This means that in the 
case the quark is produced by a hard scattering in a medium L is the length of matter measured 
from the production point. 

The subtraction of the value of the integrals at k = in fl3~i~l) eliminates the medium- 
independent factorisation contribution 11 . 

It is straightforward to simplify (j3l[). Using pa = 1/A and ( p2|) to write dt = j^Xdr, and 
defining ro = N c L/2CfX, one finds 



duo dz 7r 2 A 



ujdl a s N c Re J Su (£ dT / _ _r_\ /(£ _ xZ; r) . f *(u_ x v) 



1 r°drf(U-xV,r) ^" 



T Jo (u-xvy 



(32) 



It is convenient to express the integrals in impact parameter space. Using ( |T3D and ( |17b| ) one 
finds 



Re | dr M - - f(B, r) ■ f*(B) 



dujdz 2ir 3 A I ./n V T n / 2ix 

2m r j ~c t td is}* d2B ^ k=0 



+ — r*rf<B,T 

TO Jo 

In the small-x limit (|33| ) reduces to 



(33) 



B 2 2tt 



udl 2a s N c n . , / r\ 1- U(£ 2 ) „ ?/ „ . rf 2 5 
Re! dr 1— -±—LB-f(B,r) 



dui dz 7r 2 A I J \ T o J B 2 2tt 

1 r° , B ~ 4 d 2 5 s=0 



+ — / dr=r- f(B,T, 
2rJ B 2 7 2tt 



(34) 



Comparing with (4.24) of Ref. we see that the first term on the right hand side of is 2/3 
times (4.24). This 2/3 is exactly the 2/3 factor we discussed at the end of Sec. 4. 

The second term on the right hand side of (0), due to the hard scattering in the medium, 
was missed in Ref. 0. At first glance it might seem that this term is small compared to the 
first term because of the 1/tq in front of the integral. However, the second term is enhanced by 
a 1/B 2 compared to the first term and 1/B 2 is of the order of Tq in the dominant part of the 
integral @. 
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5.2 The induced spectrum 



It is not difficult to solve fl25f ), in a logarithmic approximation for small B 2 , and to use that 
solution in (^) to obtain the induced gluon spectrum. The details of the procedure are given 
in Ref. ||, and here we emphasise the differences which occur when x is not necessarily small. 
When B 2 is small, and this will be the case for matter long enough so that many scatterings 
occur, it is convenient to write the "potential" in (|25"D as 



l-V{B{l-x)) + l-V{B) 



In terms of v(B 2 ) used in Ref. one has 



1 

N 2 



V(-Bx) 



-B 2 u(B 2 ,x). (35) 



u(B 2 ,x) 



2 1 



x 



c w/ ) ^) 



where v, and u, have only logarithmic dependence on B 2 for small B 2 . 
fo(B), from (|T6|), becomes 



km 



Til 



u(B 2 ,x) {B} . 



(36) 

In this small- B limit 
(37) 



The only change from Ref. || is B_ — > {B} and v 

itt{B} 



u. Eqs. (5.12) and (5.13) now become 



and 



f(B,r) 



f(B,r) 



2 COS 2 U1T 



exp 



27ii{B} d 



B 2 



^ 6XP 



- —muj B 2 tan ujqT 



■-mtUoB 2 tan uoqt 



(3* 



(39) 



respectively. Parameters m and u;o are defined as before but with u replacing v. That is 



m 



1 



2C f 



IKU , K 



N r 2u(l-x) 



(40) 



Eqs. (^) and (|39|) are two useful forms of the solution to (|25|) in the approximation of treating 
ii as a constant, an approximation which should be good for small B 2 . 

If one substitutes ( |39D for f(B_,r) in the first term on the right hand side of (0), and (|38|) 
for the second term on the right hand side of (|33D , then 



ujdl 



a s N c 



du dz 7iXr 



x 2 \ „ r r° dr 

x + — Re < — 



tan UqT 



1 



^ r 



sin UqT cos Wo''" 



(41) 



emerges. The two terms in the integrand in fl4T|) correspond exactly to the terms on the right 
hand side of (|33|) . In the small-x limit the first term on the right hand side of (|4lD is smaller by 
a factor 1/3 than (5.15) of Ref. 0, this factor of 1/3 being part of the 2/3 found in Sec. 4, with 
the 2 in the 2/3 going into changing a v to a u. 
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The r-integral in ( fil]) is easily done to give 
uidl 2a s Cp 



du dz 



or 



7lL 

uidl 
duo dz 



x 2 


Mn 


sin ujot 


-In 


tan ujqTq 




Y _ 











2a s C F 



x" 

+ Y 



In IcostuoTol 



(42a) 



(42b) 



We remind the reader that 



corresponds to the quark being produced in the medium. For 



quark approaching the medium from outside the spectrum is given by the first term in the right 
hand side of (pa]). 

Now we are in a position to generalise the gluon energy spectra fl42|) for the case when the 
projectile is not a spin- 1/2 fermion but a vector or a scalar object. To this end we note that the 
basic gluon emission vertex ([I]) becomes diagonal in the gluon helicity basis, s = ±1, 



0)* 



+ 8- 



X 



1 for s = +1 
1 — x for s = — 1 



(43) 



' 1 - x)\ s - rJp \x\ r - r '\ Jp 



(44) 



In general, the structure of the gluon radiation vertex is © 

2 

x(l — x) 

where Jp is spin, and r, r' = ±1 helicity states of the projectile before and after gluon emission. 
For the quark case we have Jp = | while helicity conserves, r = r' . Taken together with ((43|), 
this brings us back to ([[]): 

u rl (p - k) 7 • e (s) * u r (p) = [ (1 - x) 1 ' 2 ■ 6 s , r + (1 - x f 2 ■ 5 s ^ r } . 

x 1 — X 



Applying ([0]) we obtain the expressions in the square brackets in (f44|) for the amplitude of gluon 
emission off a scalar particle, Jp = 0, 



[ (1 - x) ■ (<Vfi + Vi) ] , 
and for the vector projectile (gluon), Jp = 1, 

[ (l • <5 r;S + (1 — x) 2 ■ S r ^ s ) 5 r y + x 2 ■ 5 rtS 8 r - r i ] 



(45a) 
(45b) 



The last term in ( 45 bp corresponds to helicity flip of the incoming hard gluon. (The full answer 
is symmetric with respect to two outgoing gluons, x <->• (1—x), s <-> r', as it should be.) 

Adding together the squared helicity amplitudes we finally obtain the x-dependent factors 
Xj p in the gluon energy spectrum for different projectiles: 



Xx(x) 

2 

X (x) 
Xx(x) 



1 + (1 



2 

[1-x] = 

1 + (1 - x) A + x A 
2(1 -x) 



x 1 + (1 — x) 
' 2 x 

| -2(1- a;), 



x 
2 



x 



+ - — - +x(l -x)\ . 

X X 



(46a) 
(46b) 
(46c) 
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Factors fl46|) are identical in the soft limit, Xj p (0) = 1. They are proportional to the correspond- 
ing parton splitting functions. 

If the projectile corresponds to colour representation, R, which is different from the funda- 
mental representation, the colour factor Cf in the gluon energy spectrum ( fj2|) should be replaced 
by the appropriate quadratic Casimir operator, Cr. The colour factor 1/N 2 in the third term in 
the "potential" ( |35l) should be replaced by 

J_ 2C R - N c 

N? N ' 



which leads to the substitution Cp — ► Cr in the expression (|3q ) for u 



5.3 The energy loss 

The energy loss per unit length, 



dE [°° udl , 



dz /n duj dz 



can be evaluated easily if the dominant values of u ~ f^ 2 L 2 / A are such that the small-x approx- 
imation can be used when doing the u;-integral of (0). In this case 

dE a s N c fi 2 L _ x 

-!; = - — r v{T ° h (47) 



where we have used u ~ 2v in the small- a; limit. We remind the reader that in (I7[) L is the 
length of material traversed by the quark beyond its production point. 

The relative contributions of the first and second terms on the right hand side of fljlp to p7| ) 
are 1/3 and 2/3 respectively. Thus for a quark approaching the medium from outside (f|7]) is 
replaced by 

dE a s N c fi 2 L _ x 

- d7 = 12" v{T ° } ' (48) 

Comparing with Ref. [pi in the small-x limit we note that the first term on the right hand side 
of ( |42a| ) is a factor 1/3 times the expression given in (5.16) of Ref. with, of course, v replaced 
by u in the definition of Uq in the present result. 

Comparing energy loss with jet broadening @, for example for a jet produced in matter, we 
find 

_dE _ a^Nc 

dz ~ 4 P±w > (49j 

the coefficient being a factor of 2 larger than that quoted in Ref. [|J . 



6 Comparing to the method of Zakharov 

Recently B. Zakharov 611 has proposed a concise and elegant formulation for describing and 
calculating the energy loss of high energy partons in hot and cold matter. At first sight Zakharov's 
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formalism appears very different from that of BDMPS. In Zakharov's picture radiative energy 
loss of a high energy quark is described by the interaction of a high energy colour-neutral quark- 
ant iquark-gluon system with the QCD medium it passes through. No trace of a quark radiating 
a gluon as it passes through a medium and carries out multiple scatterings with that medium 
remains visible. Nevertheless, as we shall see below the two formalisms are in fact exactly 
equivalent. We begin by casting ([33|) in a form given by Zakharov. Then we shall attempt to 
explain how one can intuitively see the equivalence of the two formalisms. 



6.1 Quantitative equivalence of the two formalisms 



In this section we shall show that fl33|) can be expressed in a form identical to Eq. 4 of Ref. 
It is convenient to go back to the unsealed variables used in (|3lD, and also to use 



f(B 2 ,t 2 ) = J d 2 B 1 G(B 2 ,t 2 ;B_ l ,t 1 )f (B 1 ), (50) 

where G is exactly as given in (5.6) of Ref. || with, of course, the replacement of v by u being 
understood. Then fl33|), or (|3lT), takes the form 



G(B 2 , z 2 ; B,, Zl ) 1 \ 1; > -2-p( Zl ) a + / dz 2 p(z 2 )a fi(B 2 )G(B 2 , z 2 ; B v 0) 

where we have allowed the integrations over z\ and z 2 to go to oo. It is understood, however, 
that p(z) vanishes outside the interval < z < L. We have used z rather than i as a variable 
to emphasise the spatial dependence of the density p. In what follows it is not important that 
p be uniform in the region < z < L, as we are only going to the equations for G and not the 
explicit solution. The Green function G obeys the same differential equation as given for / in 
(|5|). Using (^) and (|35|) one finds 

and 

Using (|52b| ) we substitute the combination of dG/dzi and V 2 Bl G for the \B\u ^F- p{z\)a G in 
the first term on the right hand side of (|51|). The lower limit of the ^-integral of dG/dzi exactly 
cancels the second term on the right hand side of (|51~1) while its upper limit, at z 2 , contains 

G{B 2l z 2 ;B l)Z2 ) = 5 2 {B 2 -B X ), 

comes out ^-independent and therefore cancels due to the k versus k = subtraction indicated 
in (0). 
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Thus, we arrive at 



uodl Oi s Cp ( kN ( 



<Ijj (I:. I \2XC F 

POO 

dzi \ dz 2 p(z 2 )a 



Re ^ / d 2 B 2 / d 2 B 1 



Nc_ 
2C f 



■f^B 2 )V 2 Bl G(R z ,z 2 ;B l ,z 1 ) 



k=0 



(53) 



Now, using flSTp and ( p2a|) gives 
uidl 2a s Cp f kN ( 



duj dz 



Lit 2 \2XC f 



Ret I I d 2 B l d 2 B 2 { - l] ' { - 2] 



B\ B 2 



V 2 



d 



dzi I dz 2 \—G(B 2 ,z 2 ;B 1 ,z 1 ) + i 



2 

\2\Cf 



^im^^b^z,: 



k=0 



(54) 



The ^-integral of the dG/dz 2 term vanishes: at the upper limit, ,22 = 00 because G(z 2 = 00) = 0, 
and at the lower limit due to ^-subtraction, as before. The remaining term in (fjj) can be 
simplified by writing 



f = i{V B }.n^. 
After integrating the {Vb} and V 2 B terms by parts one finds 



(55) 



udl 2a s Cp ( RN t 



duo dz Lit 2 \2XCf 



Vl^lnS? 



Re / / d 2 B 1 d 2 B 2 / dz 2 / dz x 



1 



V| 2 -lnS 2 2 



(56) 



Using 



and 

one finds 

cud/ 2a s Cp 
(ia; dz L 

Noting that 



{V Bl } • {VbJ* = ( 1 - X + y ) V Bl • V B2 



V|-ln£ 2 = 2ir6 2 (B) 



[A-Ax + 2x 2 ] 



2XCf 



OO PZ2 



Re / dz 2 / dzi [ V Bl Vij 2 G (5 2 , z 2 ; B_ x , z x ) 
Jo Jo 



J_Bi=S 2 =0 



(57) 



(5? 



K=0 



(59) 



2AC F 2x(l - x)p 
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and that /x normalises our "impact parameter" B to physical transverse coordinate, x, so that 
f 1 Vbi = V £ , /x 2 G(B 2 , 0^,0) = ^fe-xj, 

we see that fl59| ) takes exactly the same form as Eq. 4 in Ref. [H]. 

It is not immediately clear that the subtraction terms, done at k = in our case and at zero 
matter density in Ref. ||, are identical. There is, however, a physical argument that shows that 
they should be the same. The k — >0 limit is equivalent to the u-^oc limit. However, at large uj 
the gluon is surely radiated outside the medium and since the high energy quark is produced in 
the medium, the gluon has no knowledge of the medium whatsoever. Thus, subtracting out the 
zero-density calculation is the same as subtracting out the large-tu gluons. 

In closing we note H that the two terms on the right hand side of ( flT| ) correspond to in- 
tegrations when z\ and z 2 are in the medium and when z\ is inside the medium while z 2 is 
outside, respectively. Thus a formula like (pT5|) but with z\ and z 2 restricted to lie in the medium 
reproduces the induced radiation off a high energy quark approaching the medium from outside. 
In this case the subtraction at k = subtracts out the so-called factorisation contribution fl S . 
It does not appear that this subtraction can be done in terms of a zero density limit as the 
subtraction term has a (weak) matter dependence. 



6.2 Why the two formalisms are equivalent 

In this section we describe, very qualitatively, how the two formalisms are related. We do this in 
the context of a hard scattering producing a high energy quark jet in a finite-size medium. The 
quark then radiates a gluon either in the medium of after it has left the medium. In the BDMPS 
approach one calculates the gluon emission amplitude and evolves it in time up to the time that 
the gluon is also emitted in the complex conjugate amplitude, at which time it is certain that 
the gluon will be produced and is not just a virtual fluctuation. After the gluon is emitted in the 
complex conjugate amplitude the gluon emission spectrum is determined and it is not necessary 
to follow the system any further in time. 

Suppose the gluon is emitted at t = t± in the amplitude and at t 2 in the complex conjugate 
amplitude. In evolving the amplitude, and the complex conjugate amplitude, between t\ and 
t% the quantum mechanical phase depends on the energy of the quark and of the gluon in the 
amplitude and on the quark in the complex conjugate amplitude as indicated in (ITS]) where 



the phase contribution from the complex conjugate amplitude comes with a sign opposite to 
that of the amplitude. Thus, formally, one may insert the phase from the complex conjugate 
amplitude into the amplitude by introducing a negative kinetic term in the effective 2-dimensional 
Lagrangian, and that is what is done in Refs. || |(| and [TJ] . 

As the quark- gluon system in the amplitude, and the quark in the complex conjugate am- 
plitude, evolve between t\ and t 2 there may be inelastic collisions with the medium involving 
both the amplitude and the complex conjugate amplitude. If the elastic scattering amplitude 
of a quark with a particle in the medium is purely imaginary then the total inelastic scattering 
cross section is given by the forward elastic amplitude. Then the inelastic contribution can be 
taken into account, solely in the amplitude, if one brings the quark from the complex conjugate 
amplitude to the amplitude as an antiquark i'H. 

Thus one can consider the evolution of a quark- ant iquark-gluon system in the amplitude, and 
with the antiquark having a negative kinetic energy term, as equivalent to the evolution of a 
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quark-gluon system in the amplitude and a quark in the complex conjugate amplitude. Forward 
elastic scatterings of the qqg system with particles in the medium are equivalent to the various 
elastic and inelastic reactions of the original problem. These forward elastic scatterings may be 
viewed as a two-body imaginary potential between the various pairs of the three-body system. 
This is also apparent in fl23|) and (p5|). 

Thus, formally, one can represent gluon emission in terms of the evolution of the amplitude 
for the three-body quark-antiquark-gluon system in a medium. There is a bit of a mystery in 
both approaches, and that concerns the number of independent transverse variables required to 
describe the process. One might expect there to be two independent impact parameters necessary 
to describe the three-body evolution. In the BDMPS approach one also would, in general, expect 
two impact parameters to be necessary to describe the amplitude f{U_iY--> T ) since there would 
seem to be two independent momenta, U_ and V_. In the BDMPS approach we have seen that U_ 
and V_ only appear in the combination U_ — xV_, as is apparent from (|7]), (^) and (0), so that 
only one coordinate is required. 

The structure of the potential in ( p5[ ) can be visualised as pairwise interactions between a 
quark in the amplitude put at transverse coordinate 0, the gluon at coordinate B_ and the quark 
in the complex conjugate amplitude at xB_. This same result was earlier found by Zakharov 
in Ref. ||. If x g , x_g and x g are the transverse coordinates of the quark, the "antiquark" and 
the gluon then the two-body interactions only depend on relative distances, so that the total 
momentum V is conserved: 

(1 — x) ■ Xq + x ■ x g — x_g = V — const . 

The factors (1 — x), x and (—1) correspond to the relative "masses" of the quark, gluon and 
antiquark, respectively. With the boundary conditions x g — x^ = = x — Xg at t = ti, t% one 
gets 

(1 - x) (x^ - Xg) + x (x g - Xq) = . 
Choosing x^ = one finds Xg = xXg just as in the BDMPS approach. 
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Figure 2. 




